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Abstract 

The variational principle (VP) has been used to capture the metastable states 

of a glass-forming molecular system without quenched disorder. It has been 

shown that VP naturally leads to a self-consistent random field Ginzburg- 

Landau model (RFGLM). In the framework of one-step replica symmetry 

breaking (1-RSB) the general solution of RFGLM is discussed in the vicinity 

of the spinodal temperature Ta in terms of "hidden" formfactors ^(k), 50 (k) 

and A(k). The self-generated disorder spontaneously arises. It is argued that 

at r < Ta the activated dynamics is dominant. 
PACS: 05.20.-y, 61.43-j, 61.43.Fs 
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It is a common knowledge that the transition from a hquid to a glass is mainly determined 
by extremely slow dynamic processes. On one hand, the mode coupling theory (MCT) 
predicts a pure kinetic transition temperature, Tc , where the density fluctuations are frozen 
and the system becomes nonergodic |T[]. On the other hand, there is a strong trend in favour 
of the idea that the appearance of metastable states in the free energy landscape significantly 
determines the dynamical arrest at some transition temperature Ta [H^- This consideration 
was based on both Potts and p- spin glass models, in which the quenched disorder 
was given. The authors drew the conclusion that these spin-glass models without reflection 
symmetry are also relevant for the structural glasses, where the disorder in a sense should 
not be quenched but self-generated 0. 

Recently a great deal of work was done to incorporate the ideas and notions of the mean- 
field theory of the spin-glasses into the theory of the structural glasses [[7|-p!l[]. The main 
underlying idea is related to the free-energy landscape. In order to quantify this landscape 
paradigm several copies ( or "real replicas" ) of the same system have been introduced. In 
doing so a reference replica ?/ is a typical configuration distributed with a Boltzmann - Gibbs 
measure, whereas a replica x is affected by the influence of the replica y through some short 
ranged attractive potential. It is assumed that for a fixed value of y, which plays the role 
of the quenched variables, the value of x thermalized at the temperature of the experiment. 
The small positive coupling constant e is imposed which forces the variables x to stay near 
the y variables. If this correlation between x and y variables persists even at e — > the 
system is in a glassy phase. As a main result the expression for the "effective potential" has 
been calculated which shows a nontrivial solution for the "overlap" q{x,y) [^j-j^. 

In this letter we shall present a different general method in order to treat the metastable 
states in an arbitrary glass-forming molecular system without the quenched disorder. For 
this purpose we reduce the initial problem to a self-consistent random field Ginzburg-Landau 
model (RFGLM) by making use of the Feynman variational principle (VP). We shall show 
that the VP automatically yields the standard replica structure in phase space with the final 
limit n — for the number of replicas. This RFGLM is investigated by implementation of 
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more or less standard field-theoretical tools [p!^ -[T5[]. The replica symmetry breaking (RSB) 
manifest itself as a 1*^* order phase transition in the framework of a one-step RSB (1-RSB) 
scenario and the self-generated disorder shows up. Some comparison with the "real replica" 
method [[7|-p!l[| will be given below. 

We shall start with a low-molecular system of particles with coordinates r^^^^ , (p = 
1, N) via a pair potential V^(r). The partition function at temperature T = 1/(3 is given 
by 



N ( n N 



Z= ff[ rf'r(P) exp i Yl ^(^^""^ - ^^™^) [ • (1) 

p=l I, p,m=l ) 

A physically meaningful representation of eq. d^) in terms of the collective density, p(r) = 
Ylp=iK^ - r^^^), is given by 

Z = j Z}p(l)exp|iy[p]-| j rflrf2p(l)\/(l-2)p(2)|, (2) 

where the entropy of the system is related to the Jacobian of the corresponding transforma- 
tion 



N 

W[p] = log / JJdl^*')^ 



N 



N 

p=l 



p=l 

and the short - hand notation r^^^ = l^^^ has been used. An alternative appropriate repre- 
sentation of eq. (0) can be obtained in terms of a field tpi^^) conjugated to p(r) which appears 
as a result of the functional Fourier transformation 

exp{K[^]} = J Z}p(l)exp|iy[p]-z J rflp(l)7A(l)| , (4) 

where K[ip] is the cumulant generating functional for the free system. The functional ex- 
pansion for K[iIj] can be simply obtained in an explicit form. The invertion of eq. @ and 
substitution in eq. (|I]) gives 

Z = I D^exp|i^[^]-i j rflrf2V^(l)[/?V]-i(l-2)^(2)|. (5) 

As a result, eqs. (|^) and provide two equivalent representations of the partition function. 



Our main objective in the following is to calculate the branch of the free energy which 
appears in a supercooled (below the melting point) liquid and which corresponds to the 
metastable states. To this end we recall that in a glassy phase only the component average 
free energy 

F = + TE (6) 
is physically meanigful |T^. In eq. (^ the canonical free energy = —TlogZ and the 



function S > is called the complexity |T6[. On the other hand it is known that the 



Feynman VP, when implemented for the partition function (|^), yields an upper limit for Fc 
, i.e.Fc < Fvp. As the basic assumption we take therefore 

F = rmni^F,+ (^^ j dld2p{l)V{l - 2)p(2) - W[p] - S[p\J | , (7) 

where the trial (Gaussian) Hamiltonian 

S[p] = \ I dld2{G-\l,2)p{l)p{2) -2G-\l,2)p{l){p{2)),} . (8) 

Here the expectation value (p(l))s and the correlator G(l, 2) = {6p{l)6p{2))s are considered 
as independent variational functions. The assumption (|^) means that the variational free 
energy Fyp takes its minimal value on the ensemble of components (or pure states). It will 
be justified a posteriori by the appearance of the replica structure. The minimization in 
eq. (1^ with respect to (p(l))s and G{1,2) leads to 

(p(l)), = J d2[pV]-\l,2) (^W^h) , (9) 
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1,2) , (10) 



where the "quenched"?/'- expectation value reads 



/ 5 \ _ JD^P^mew{m+^Idlpilm)} ^ 

\6p{2) ^^V, \ jDijexp{m+zjdip{mi)} / ' ^ ^ 

and the "quenched" ip - correlator 



/ / fD^P5m5i;{2)exp{K[i;]+iJdlp{mi)] ^^^^ 



\6pil)6pi2) ^^V^ \ jD^exp{K[^]+ifdlp{l)ij{l)} 

arises naturally in this representation. As a result we have reduced the problem to a self- 
consistent RFGLM . We have used the term "quenched" because the Gaussian field 
plays the role of an "external quenched" field in eqs.(]TT|)- Its moments however should 
eventually be determinated self- consist ant ly from eqs. (^-(|IU|). 

One should not be surprised that Nature reveals the replica structure through VP. It is 
well known ( see e.g. Sec. 4 in that the rigorous treatment of statistical thermodynamics 
yields only thermal equilibrium states. In order to capture the metastable states in the 
framework of statistical thermodynamics one has to constrain properly the phase space of 
the considered system. The most famous example is the van der Waals - Maxwell metastable 
loop, which exists in the mean - field approximation and becomes fiatter if the clusters 
formation are taken into account [|T^. It is the VP in our case which implements this 
constraint. 



The "quenched" moments (|Tl|)-(|T2D can be calculated by the standard way ||^ from the 
replicated partition function 



(Z")JG(l,2),(p(l))J 



J]DV>,exp|-- / dld2j2 [po5afc5(l,2) + G(l,2)]^„(l)^fe(2) 



a=l a,b=l 



/n n 
rfl(p(l)).5^^„(l) + 5^Aq^j} , (13) 
a=l a=l 



where po is an average density and -^[V^a] is the anharmonic part of -A['?/'a]. As has been 
shown in ref. |T^-|TB[ the effective Hamiltonian such as in eq. (|TB|) may lead to the replica 



symmetry breaking (RSB) which corresponds to a structural glass transition. 

It is important to mention that for potentials V{r) with the infinite range interactions we 
can directly expand the effective Hamiltonian in eq. (0) around the saddle point solution p(l) 
up to the second order. The calculations show that this next to the mean-field approximation 
and the VP merge and that both become exact, i.e. F = F^. In this case the glassy phase 
does not appear. This conclusion for the present realistic model (|l|) is in agreement with the 



result for a 0{M) - model with a quartic interaction in the large M limit as well as for 
the particles on a hypersphere [0. For the model given by eq. (|I|) with the infinite range 



of interaction the self-generated disorder is not generic and only crystallization can be 
expected. It can be shown also that the Langevin dynamics of the pure ( i. e. without a 
quenched disorder ) system with the screened Coulomb potential V{r) = {^/N) exp(— 7r)/r 
( where the interaction parameter /i > , the inverse screening length 7 oc N~^^'^ and 
the number of particles >> 1 ) does not reveal any glass transition. Despite ongoing 
discussions [^] the glassy dynamic behavior of the pure models is still debated. 



Now we shall mention briefly the technical tools to calculate the free energy 
•>^{^afe(l, 2), as a functional of the expectation value {ipai'^)) and the correlator 

^a6(l,2) for the rephcated field theory given by eq. (|1^). To this end, as in ref. [|l^,|l^, we 



have used the second Legendre transformation. The detailed analysis shows that the Hartree- 
like contribution into the so-called generating functional of all 2-irreducible diagrams is zero. 
The first nontrivial contribution comes from the second-order terms in powers of vertices in 
eq. (P^. The resulting free-energy functional has been parametrized in the framework of 
the Parisi Anzatz 0: ^^6(1, 2) — > {g{k), g(k,x)} where the formfactor ^(k) corresponds 
to the diagonal and (^(k, x) to the off-diagonal elements of the hierarchical matrix Qabi^,'^)- 
The wave vector k- dependence is retained in the spirit of ref. ||2^ . In one can also find 



the beneficial formulas regarding the algebra of the hierarchical matrices. In the off-diagonal 
formfactor g(\<.,x) the argument x G [0, 1]. 

We have restricted ourselves to the simplest case of 1-RSB scenario p . This means that 
the off-diagonal formfactor (yf(k, x) consists only of two pieces: (^(k, x) = go(k) at x < Xc and 
g(k,x) = gi(k) aX x > Xc , where Xc is the break point. It has been shown before that for 
both Potts and p-spin glass models P,|23[] the 1-RSB scenario is generic. 

It is natural to define the gap A(k) = gi(k) — (7o(k) as a non-ergodicity (order) parameter 
of the required glassy phase transition. As a consequence of these calculations the free-energy 
functional {^(k), (yfo(k), A(k)} has been obtained. The extremization of it with respect to 
^(k) and go(X) at A(k) = yields the equations for the replica symmetric (RS) case 
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1 1 I ^^o(k) , ^ 

+ r~/1 N /I M9 + Po 



^(k)-c(k) m-9o{^) m)-9om' 
+|y"^(k-ko^(ko-| 1 ^(k-ki-k2MkiMk2) = o 



ki ki,k2 

5(o(k- ki)5(o(ki 



^o(k) , 1 pI 



[^(k)-^?o(kP m-C{k) 2 

ki 

y (7o(k-ki-k2)^7o(ki)(7o(k2) = , (14) 



ki,k2 



where J = f df'\<i/{2TiY and we used the approximation —(3V{\<i) = C(k) (C(k) is the 

k 

direct correlation function ||2^). For simphcity we have consided the homogeneous case, 
(p(l))s ~ const, to ehminate the hnear term with respect to ip in eq. (p!3|). 

In order to grasp the glassy (or 1-RSB) solution , {g, qq, A} has been expanded up to 
the 4-th order with respect to A(k) around the RS-solution (|T^. As a result the increment 
of the free energy connected with the non-zero order parameter A(k) is determined by the 
Landau expansion 

J'i{A} = - J r2(qi,q2;x,)A(qi)A(q2)- j (qi, q2, qs; A(qi) A(q2) A(q3) 
qi,q2 qi,q2,q3 

- J r4(qi,q2,q3,q4;a;c)A(qi)A(q2)A(q3)A(q4) + ... (15) 



qi>q2,q3,q4 



The form of the expansion coefficients is determined by the RS-solution of eq. (113) as it 



will be given in an extended paper I^Sj. Since is generally speaking nonzero the transition 
is discontinuous. As in ref. our free energy is proportional to the self-overlapping 

fraction : jFi oc — (1 — Xc)- In this has been argued in the metastable states 

appear first at the spinodal point Ta which is determined by 

,„„^/^l .0 . (16) 



xv— >1 



A=A 



5A(q) \l-xj 
provided that the eigenvalues of the matrix 



(17) 

A=A 



are nonnegative. 

It is an open question whether the MCT transition temperature is related to T4 . For 
spin models with quenched disorder has been shown that these temperatures coincide . 
In any case at T < the dynamics is dominated by the activated process PJ26|. 

It is interesting that within the "hidden" correlators ^(k), (70 (k) and gi(k) only ^(k) 
is related to the density correlator: G{k) = [g(k) — C(k)]^^ (see eq. (p!Ol)), so that the 
discontinuous transition does not shows up through the density correlator G'(k) which is 
measured by an experiment. Nevertheless one can expect that below the spinodal point T4 
the activated dynamics becomes dominant and manifests itself in the time-dependent density 
correlator G{]s.,t). We shall return to this question in the near future. The equations for 
the RS-case ([T^ ) should be solved numerically first and then one can solve eq.(|l^)- {^i7\ ) for 
the order parameter A(k). 

The glassy state appears first as a metastable one (— J^i {A} > 0) but upon further 
cooling it can become stable: —J-'i {A} < (we recall that the conditions (|16|)- (0) deter- 
mine the local maximum as usual in spin - glasses |^). This possibility has been discussed 
in ref. where on a basis of the Ramakrishnan - Yussouff free energy functional many 



inhomogeneous ( glassy ) density configurations have been found. On the contrary, in spin 
- glasses the glassy state is thermodynamically metastable for Tk < T < Ta { Tk is the 
temperature where the complexity S becomes non-extensive) p|-^,^|. 



We would shortly like to compare our method with the mentioned "real replica" ap- 
proaches |7HTT|. The main difference is that in this method we deal with the conventional 
replicas where the "zero replica limit" and RSB are already fairly standard 0. This leads, 
among other things, to a close link with the generalized spin-glass models results P-^,p7 



Our consideration based on the well defined replicated field theory (see eq. (|1^) ) which 
needs as "an input" only information about the direct correlation function . Conversely, 
in the "real replica" approach one has to choose one or other simple liquid models ||2^ 
(e.g. the replicated hypernetted chain (HNC) approximation in ref. [^[1^1 ). Finally, our 
considerations can be directly generalized to the Langevin dynamics by making use of the 



MSR- generating functional method p0|j28[] , whereas it is rather unclear what a dynamical 
counterpart of HNC in the simple liquid theory is. 

In summary , we have suggested a general method to deal with the equilibrium structural 
glass transition. Making use of the VP naturally leads to the self-consistent RFGLM and the 
RSB - transition at some spinodal temperature T^. In the framework of 1-RSB scenario the 
equations for the set of "hidden" formfactors ^(k), go(k) and A(k) are derived. Although 
the static density correlator G(k) does not show the RSB its dynamical counterpart could 
be considerably affected by this transition through the onset of the activated dynamics. 

The authors have benefited from discussions with J.Baschnagel, K. Binder, W. Kob, 
A.Heuer , K.Miiller-Nedebock and gratefully acknowledge financial support by the Sonder- 
forschungsbereich 262. 
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